Experimental Verifications of the Casimir Attractive Force 

between Solid Bodies 



S.K. Lamoreaux 

University of California, Los Alamos National Laboratory, Physics Division 23, M.S. H803, Los 

Alamos, NM 87545 
(February 1, 2008) 

(Submitted for the Casimir 90th Birthday Issue of Comments on Atomic and Molecular 
Physics) 

Abstract 



(N 



> 

. A brief review of the recent experimental verifications of the Casimir force 

C . between extended bodies is presented. With modern techniques, it now ap- 
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pears feasible to test the force law with precision better than 1%; I will address 
the issues relating to the theoretical interpretation of experiments at this level 
of accuracy. 
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I. INTRODUCTION 

The force between uncharged conducting surfaces, the so-called "Casimir force," has 
been described as one of the least intuitive consequences of quantum electrodynamics. For 
conducting parallel flat plates separated by distance r, this force per unit area A has the 
magnitude Q: 

7r 2 fir 1 

F{d)/A = = 0.013-dyn(/.m)7cm 2 . (1) 

This relationship can be derived by considering the electromagnetic mode structure between 
the two plates, as compared to free space, and by assigning a zero-point energy of ^hto to 
each electromagnetic mode (photon). The change in total energy density between the plates, 
as compared to free space, as a function of separation d, leads to the force of attraction. 

The only fundamental constants that enter Eq. (1) are h and c; the electron charge e is 
absent, implying that the electromagnetic field is not coupling to matter. The role of c is to 
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convert the electromagnetic mode wavelength to a frequency, while % converts the frequency 
to an energy. 

The term "Casimir effect" is applied to a number of long-range interactions, such as those 
between atoms or molecules (retarded van der Waals interaction), an atom and a material 
surface (Casimir- Polder interaction), and the attraction between bulk material bodies. The 
latter effect is generally referred to as "the" Casimir force and depends only on the bulk 
properties of the bodies under consideration; I will limit my discussion to the latter effect. 

For real materials, Eq. (1) must break down when the separation d is so small that 
the mode frequencies are higher than the plasma frequency (for a metal) or higher than 
the absorption resonances (for a dielectric) of the material used to make the plates; for 
sufficiently small separation, the force of attraction varies as l/rf 3 , as discussed in particular 
by Lifshitz [f2| . In analogy with the attractive forces between atoms, the force in this range 
is sometimes referred to as the London- van der Waals attraction, while the l/d A range is 
referred to as the retarded van der Waals (Casimir) interaction. For the Casimir effect, the 
crossover distance between the regimes is d ~ 100 nm, much larger than atomic spacings 
in the materials, so it still makes sense to describe the materials by their bulk properties 
(index of refraction); the 1/d 3 vs. l/d A interaction is in this case due to the truncation of the 
mode frequencies that are affected by the changing plate separation. Therefore the crossover 
between the two regimes appears to be of physically different origin than in the case of the 
attractive forces between isolated atoms. 

The Casimir effect and its calculation represent an electromagnetic waveguide problem, 
where imperfect materials are used in the construction of the waveguide. The zero-point 
fields are those associated with the waveguide; these modes do not exist in free space, so the 
idea that the Casimir effect represents a negative energy density compared to free space is 
incorrect. 

Given that the distances where the force of attraction is sufficiently strong to be ex- 
perimentally detected are d ~ 1000 nm or less, an accurate theoretical description of an 
experimental system must take into account the real material properties as will be discussed 
in Sec. Ill below. 

Finally, it is interesting to note that there have been only a dozen or so published 
experimental measurements of the Casimir force, to be compared with the hundreds and 
hundreds of theoretical papers on the subject. Perhaps very few doubt the strict validity of 
Eq. (1) or its modification for real materials ||. Because of the unavoidable uncertainties 
in bulk material and surface properties, verification of Eq. (1) as a test of QED will likely 
always be inferior to measurements of the Lamb shift or g — 2 of the electron. 

II. EXPERIMENTS 
A. Overview 

The experimental situation as of 1989 has been reviewed by Sparnaay in the volume 
prepared in the honor of Dr. Casimir's 80th birthday || . Since then, two experiments have 
been performed, both with significantly better accuracy than had been previously obtained. 
These two experiments were based on a torsion pendulum balance Q and on atomic force 
microscopy (AFM) ||. 
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The recent experiments employed techniques that were developed in particular by van 
Blokland and Overbeek || in the measurement of the attractive forces between metallic 
films. Measurements between metallic films pose difficult problems as compared to dielectric 
films for which optical techniques can be used for alignment and distance measurements. 
In the case of metallic films, the distance is determined by measurement of the capacitance 
between the plates. Alignment is simplified by making one plate convex, in which case the 
geometry is fully determined by the radius of curvature at the point of closest approach, 
and the distance between the plates at that point. This technique was first put forward by 
Deryagin |7| and has found broader application as the Proximity Force Theorem ||. For 
this geometry, Eq. (1) becomes 

7r 2 1 he 

F(d) = 2vRE(d) = 27T— (2) 
v ' K ' 240 3d 3 V ' 

where R is the radius of curvature and E(d) is the energy per unit area that leads to the 
force in Eq. (1). It should be noted that the plate area does not enter into Eq. (2). The 
two recent experiments each used one convex and one flat plate. 



B. Torsion Pendulum Experiment 

Because I was directly involved in this experiment, I will give an anecdotal account of 
it; specific details, schematics, and data can be found in [|J. This experiment was started 
in 1994 as an undergraduate project. Dev Sen, a junior at the University of Washington at 
that time, was working with our atomic physics group; he heard of the Casimir force in his 
Electromagnetism course and asked if we might be able to measure it. With no knowledge 
of previous experimental work, with the exception of the Sparnaay paper of 1958 |J (which 
is routinely quoted as "the" verification of Eq. (1)), I proposed an electrostatically balanced 
hanging (Cavendish) torsion pendulum technique; this is different from the horizontal torsion 
balance used by van Silfhout [|1(J (which, at the time, we did not know about) particularly 
because the angular force constant a is at least two orders of magnitude smaller for the 
hanging pendulum. For short times (a few seconds of averaging) our principal noise source 
was thermal fluctuations (e.g., Brownian motion); in this limit, the signal to noise scales 
as \j\fa. It appears that this noise was a major limitation of our experiment, particularly 
when measuring the force at large plate separation. 

We constructed the apparatus with junk found around the lab, and the total material 
construction cost of the apparatus was around $300. The torsion pendulum fiber was about 
60 cm in length; mounted on one arm of the pendulum was a flat plate coated with a Cu/Au 
film, while the other arm served as one electrode of a differential capacitor system that was 
used in a feedback system to keep the torsion pendulum at a fixed angle. The second plate 
used for the Casimir force measurement was mounted on a precision mechanical mount and 
could be moved about 10 /xm by use of piezoelectric transducers. 

In our first attempts to measure the Casimir force, two flat plates were used. We were 
never able to properly align the plates, and we eventually gave up when Dev graduated. 
Some months before I moved to Los Alamos, I decided to give the apparatus one last try. I 
had the idea to replace one flat plate with a convex plate; I also improved the vacuum and 
electronic circuitry. 
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In June of 1996 the apparatus became operational; when the two plates were moved 
into close proximity, the increase in force over the residual electrostatic interaction was 
easily evident. Unfortunately, the apparatus was not stable enough to make a measurement; 
about 20 minutes after adjusting the distance between the plates so that the separation at 
closest approach would be of order 1 /im or less, the apparatus would drift until the closest 
approach would increase to over 5 /im. I eventually discovered the source of the drift: the 
concrete floor of the Physics Department basement would distort while I was standing near 
the apparatus during adjustment, and the floor would relax to its unstressed position over 
a 20 minute period. This effect was dealt with by standing in a different position, rather far 
from the apparatus, while making the initial adjustments; the experiment thereafter became 
physically (in an anatomical sense) painful. 

Data was taken over a one month period, after which the apparatus was calibrated. The 
final measurement was the radius of curvature of the Cu/Au coated convex plate (lens) 
in which I made an error. To make this measurement, I used the angular deviation of a 
laser beam reflected from the lens surface as the lens was translated normally to the beam. 
Measurements near the center yielded R = 13.5 ±0.5 cm, while translations across the entire 
lens yielded R = 11.3 ± 0.1 cm; it did not occur to me that the lens might be aspheric (it 
was a high quality, precision optical lens of unknown origin or application that was coated 
with a thin layer of Au on top of a relatively thick Cu layer). I assumed the latter number 
was more accurate, and the data was described by Eq. (2) with no corrections for finite 
conductivity; I mistakenly convinced myself that the corrections should be less than 5%. At 
this point, my work was rudely interrupted by my relocating to Los Alamos, and I published 
my results in Physical Review Letters. 

Subsequently, a number of theorists expressed surprise that no correction due to finite 
conductivity was required JTTJ, and several expressed an interest in doing a proper calcula- 



tion. None did the calculation, so I figured out how to do it myself ||12|| ; I knew that the 
plasma model correction to first order was not applicable (see Sec. Ilia); being upwards of 
a 30% correction, none of the uncertainties in calibration were large enough to encompass 
it. However, the variation in the measurement of R was in the back of my mind, and in 
October of 1997 I made a trip to Seattle and retrieved the lens. Using a mechanical gauge, I 
measured R = 12.5 cm in the region where the Casimir force measurement was made. The 
accurate finite conductivity correction is small enough that its effect is indistinguishable 
from a calibration error, or a mismeasurement of R. In the end, my data did not agree 
with the assumption that the film was pure Au, but was better described by a pure Cu 
film jnj. This result is not so surprising considering that the properties of the film depend 



on preparation technique, purity, and the possibility that the Au diffused into the Cu layer 
significantly. A more accurate calculation of the finite conductivity effect would require a 
direct measurement of the complex permittivity of the films. 

In hindsight, it is remarkable that the torsion balance experiment, which was intended 
as a demonstration, worked as well as it did. The improvement over previous measure- 
ments is due to a number of factors, including the high sensitivity of the hanging torsion 
pendulum and its lack of mechanical hysteresis, larger measurement distances so vibration 
and mechanical instabilities were less important, improved piezoelectric transducers, and 
automated data collection so that large amounts of data could be analyzed and averaged. 

Much improvement over the present accuracy obtained by this technique is unlikely. The 
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apparatus was rather unwieldy with its enormous vacuum can and its susceptibility to tilt. 
The length of the torsion fiber might be significantly shortened, reducing both the intrinsic 
sensitivity (bad) and sensitivity to external perturbations (good); however we must bear in 
mind that a factor of ten improvement in sensitivity only extends the measurement distance 
by factor of about two. 



C. Atomic Force Microscopy Experiment 

In his 1989 review ||, Sparnaay discusses the possibility of using atomic force microscopy 
(AFM) to measure the Casimir force; AFM had just been invented at that time fL4| . It was 



not until late 1998 that results from an AFM Casimir experiment were reported by Mohideen 
and Roy ||. 

In this experiment, an Au/Pd + Al coated, 0.3 mm polystyrene sphere is attached to 
an AFM cantilever. A similarly coated optically polished sapphire plate was attached to a 
piezoelectric transducer and brought near the sphere. The attractive force was determined 
by reflecting a laser beam from the cantilever tip; the displacement of the laser beam on 
a pair of photodiodes produced a difference signal proportional to the cantilever bending 
angle. 

The sensitivity of the apparatus was such that the absolute force could be determined 
with a fractional error of 1% at d = 100 nm, and about 100% at 900 nm. 

The use of AFM to measure the Casimir force might be a real breakthrough; this is 
because the AFM technique is very stable and reproducible. Unfortunately, it is limited 
to measurements short distance where there are significant theoretical uncertainties in the 
interpretation of the data. For the first time testing Eq. (1) or its modification for real 
materials to better than 1% accuracy appears possible. As described in 0, it is anticipated 
that a factor of 1,000 improvement in sensitivity appears possible, which would extend the 
separation where the Casimir force can be measured to 1% accuracy to about 1 /im. At this 
distance, the theoretical uncertainties associated with the corrections for real materials, as 
described in the next Section, become much less important. 



III. CORRECTIONS 

A. Imperfect Conductivity 

Equation (1) must break down when the plate separation is so small that the mode 
frequencies being affected when d is varied are above the material resonance or plasma 
frequencies. In the case of a simple metal, the real part of the dielectric constant can be 
approximated by 

eV) = 1 " (3) 

where u p is the the plasma frequency and is proportional to the effective free electron density 
in the metal. It is convenient to introduce the plasma wavelength, X p = 2irc/uo p . Corrections 
to Eq. (1), expanded in terms of \ p /d, have been calculated to first order by Haergraves 
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15 1 and by Schwinger et al. [16fl, and to second order by Bezerra et al. JT7| For flat plates, 
the corrected force can be written in terms of a Eq. (1) with a multiplicative factor, 



F\d) = F(d) 



;| 8 X p 120 

37T d 



Xr 



4tt 2 I d 



(4) 



This equation is only valid for X p /d << 1; unfortunately, the Casimir force is large enough 
to be accurately measured experimentally only in the range X p / d ~ 1 or larger. We are also 
faced with the problem that Eq. (|3]) is only approximate. 

It is, however, possible to very accurately determine the attractive force as a function of 
plate separation by numerical calculation, provided we know its complex permittivity as a 
function of frequency: 



e Iuj 



e (u) + ie (uj) 



(5) 



where e', e" are real. With this information, the permittivity along the imaginary axis can 
then be determined by use of the Kramers- Kronig relation, 
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dx + 1. 
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This can be used in the Lifshitz expression for the attractive force 
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(7) 



where s = ^Je{ix/p) — 1 + p 2 . The numerical calculation for the attractive force between 
Au, Al, and Cu plates has been recently published fl2]| , and significant deviations from Eq. 
(|) were found. In particular, for Al with d ~ 100 nm Eqs. (§) and ([[]) differ by about 5%; 
one should note that including the third order correction to Eq. ([|) worsens the deviation. 
However, these calculations should be considered in light of the notorious variation of bulk 
and surface properties of materials due to preparation technique, purity, etc. 0,113 • 



B. Surface Roughness 

From the earliest experiments, it was realized that surface roughness would lead to an 
increase of the apparent Casimir force and therefore cause systematic errors in measurements 
aimed at verifying Eq. (1). Such effects were observed by van Blokland and Overbeek H; 



roughness has been discussed theoretically by van Bree et al. |[L9|| , and more recently in |20 

For high-quality optically polished surfaces, the roughness RMS amplitude A is usually 
of order A = 30 nm or less. For a 1/d 4 attractive force, the correction to Eq. (1) can be 
written 
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(8) 

The correction for the recent torsion balance experiment, at the point of closest approach, 
is about 1%, while for the AFM experiment, it is about 30%. 

The roughness correction was derived in the context of a 1/d 4 force law (this can be 
easily modified for the spherical plate 1/d 3 case). However, the finite conductivity correction, 
particularly as given by Eq. effectively has terms containing 1/d 5 and 1/d 6 . In principle, 
the roughness correction should be done for each power law separately, or the average force 
determined from the accurate calculation, Eq. (|?p. One should also bear in mind that the 
simple geometrical averaging procedure isn't exactly correct; a complete treatment would 
involved solving the appropriate electromagnetic rough boundary problem. However, the 
geometrical averaging is correct so long as the period of the roughness is larger than the 
separation between the plates. 



F'{d) « F(d) 



a. 



C. Effect of Thin Films on the Plate Surfaces 

Either intentionally (Au evaporated onto an Al or Cu coated substrate) or accidentally 
(formation of oxide layers) every Casimir force measurement has made use of mono- or 
multilayer coated plates. The calculation of the force for a general film configuration has 



been given by Spruch and Zhou |[21|| . The problem that I will consider here is very simple 
so I will outline its solution in some detail. 

A simple geometry that illustrates the effect of a thin material film is shown in Fig. 1; 
one of two identical perfectly conducting flat plates is coated with a thin layer (thickness a) 
of a real substance (Au, for example), and the separation between the perfectly conducting 
surfaces is d + a. This simplified problem will allow us to determine the qualitative effect of 
a thin film. 

Milonni has presented a complete calculation for the case of materials with no absorption 
and this can be easily adapted to the multilayer case (this calculational technique 



was first described in |23| and elaborated upon in ||24|| ). Following Milonni, we consider the 
case where the electromagnetic wave propagation vector in the three materials (vacuum, 
film, perfect conductor) 



K- = e-e{^ (9) 

where k is a real number, and % — 0, 1, 2 with (eo(^) — 1) representing the space between 
the plates and 2 (e2(<^) = oo) the perfect conductor, and we require Re(Ki) > 0. For the 
case considered here, we must allow for a complex e\ in which case the KiS can be complex. 
In the Appendix, the use of the techniques described by Milonni for the case of absorption 
(complex e) is justified. 

There are two type of solutions to the wave equation, one with electric vector parallel to 
the surfaces (with arbitrary orientation which we choose as the y axis), e y (z) and one with 
electric vector perpendicular to the surfaces (along the z axis), e z (z). The wave equation is 

^-e y , z (z)-Kfe y , z (z) = (10) 
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and the boundary condition for e z are (1) de z /dz and ee z are continuous, while for e y they 
are (2) de y /dz and e y are continuous (at the conducting surfaces, e y = and de z /dz = 0). 
Ignoring unphysical exponentially growing solutions, we have 

< z < d- a 
d — a < z < d 

(11) 

where the =F sets e y = or de z /dz = at the conducting boundary located at z = 0. We 
therefore have two sets of linear equations involving A, B, C for the two cases. The condition 
for non-trivial solutions of these equations is that the determinant of the coefficient matrix 
is zero, yielding the following two expressions: 

f v ( u, k, d) = 

[(e^ a + l)K l + (e^-l)K } 



e y,z\ z ) 



de z /dz 



A(e Koz =f e~ Koz ) 
Be KlZ + Ce~ KlZ 
z = d 
z = d 



f g ( u,k,d) = 

-4- ^ (,.A(^ 2K ^ a 4-1WJ 

- I- (13) 



[(1 - e 2K ^)K x + ei (cu)(e 2 ^ a + 1)AT ] 

The zeroes u> nyjnz (k, d) of f y>z determine the allowed mode "eigen" frequencies, 

The zero point energy associated with the plates is determined by assigning energy hu/2 
to each mode; 

E(d) = £ [Huj ny (k, d)/2 + hcu nz (k, d)/2\ (14) 

n,k 

(in general, the eigen-frequencies are complex, but the imaginary parts cancel as discussed 
in the Appendix), where the sum over k, in the continuum limit, becomes an integral, 

E-f^)V 2 ^ ( 15 ) 



where L is the transverse dimensions of the plate in the x, y directions. The theory of complex 
function can be used to evaluate the sum over eigen-frequencies; specifically, according to 



the argument theorem |25| , |26 



r Q^dz = N-P (16) 
2ni Jc f(z) K ' 

where C is a closed path in the complex plane, N and P are the number of zeros and poles 
within C, respectively, and the path is counterclockwise. The argument theorem can be 
modified to give the sum of the zeroes and poles: 



S 



1 



Jc f z) ^ 



2ni Jc f(z) 



Z; 



2; 



f(Zi) = QO 



(17) 



Furthermore, f'(z)/f(z) = d/ d z (log f(z)). The eigen-frequencies of physical interest lie in 
the right half plane; integrating along the imaginary axis from oo to — oo and closing the 
path with a semi-circle at infinity around the right half plane (see |22] for details), and 
integrating by parts gives 



E(d) 



HL 2 
8^2 J(l 



/oo 
di [l°g$X£, k, d) + log# z (£, k, d)] 
-oo 



where we have set u = i£, with £ real, g y , z (£, k, d) = f x ,y{^i k, d), and 

Ki = k 2 + e^O^/c 2 . 



(18) 



(19) 



Finally, we note that the poles of Eqs. ( |l2|) and (^) do not depend on d because it only 
enters in the multiplicative exponential; therefore, Eq. (pT8|) gives the zero point energy up 
to an additive constant, while the force per unit area is given by 



47T 2 JO 



k dk d£Ko 



+ 



9y(tk,d) g z (i,k,d) 



(20) 



where (possibly) non-physical ^-independent terms are omitted. 

F(d) for a Au film 35 nm thick is shown in Fig. 2; ei(z£) was determined from tabulated 



optical constants as described in 12 



IV. CONCLUSION 

In the last few years, two new experiments to verify the existence of the Casimir force 
between solid bodies have been performed and have shown a new level of accuracy; in 
particular, atomic force microscopy technique offers great promise for testing Eq. (1) and 
its modifications for the case of real materials to better than 1% precision. In particular, 
if the sensitivity of the AFM experiment can be increased substantially, the measurement 
region can be extended to larger separations where the theoretical uncertainties discussed 



above are substantially reduced. Perhaps even the effect of finite temperature |27| will be 
measurable in the not-too-distant future. 



V. APPENDIX 



As discussed by Milonni [ 22| , the contour integration technique used to sum the eigen- 
frequencies is technically correct only if the eigen-frequencies lie on the positive real axis. 
There has been much commentary on the extension of this technique to absorptive mate- 
rials (i.e., complex permittivity) in which case the eigen-frequencies are complex. Barash 
and Giszburg |24[] introduce the idea of an auxiliary system to account for the complex 
permittivity, with the fundamental eigen-frequencies real. However, the contour integration 
method gives the correct answer quite simply; the mathematics does not know about the 
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auxiliary system. I would like to offer a non-rigorous explanation for why the technique 
works in the case of absorption. 

First, for a generalized permittivity | 2"B|| , e(— uj*) = Therefore, the eigen- 



frequencies for the general boundary problem case occur in pairs, oj = ±u/ + iu" . We 
further take the case where Re(Ki) are either all positive or all negative which follows from 
continuation; as u — ► oo, all the K{ become equal because e%{oS) — > 1. In the case that 
Re(Ki) > 0, representing exponentially damped surface wave, the eigen-frequencies lie in 
the lower half plane; therefore, e~ luJt is damped exponentially in time. For Re(Ki) < 0, the 
eigen-frequencies lie in the upper half plane and represent solutions growing exponentially 
in time and space. Clearly, the contour integration method as described above should not 
work without justification regarding branch cuts etc. 

The way to see around this problem is instead of considering the eigen-frequencies, one 
can consider the corresponding eigen- wavenumbers Ki, and K in particular. The eigen- 
wavenumbers occur in complex conjugate pairs K = K' ±iKq (as can be seen from Eq. ([|) 
and the properties of e(uS) in the complex plane) and by definition the Ki are in the right 
half plane for the exponentially damped solutions. Furthermore, if we write the determinant 
function in terms of Kq, and using the fact that KqcIKq = —udu/c 2 , we find 

Jc f{uj{K ) Jc f{u{K )) 

and we can replace dK^ by du because the path is arbitrary in the complex plane (note 
that in the case of no absorption, the eigenvalues for Kq lie on the imaginary axis while 
those for uj lie on the real axis, and the contour must be adjusted accordingly). Because the 
eigenvalues for K occur in conjugate pairs, the left side of Eq. (ET) is real; therefore, the 
sum in Eq. ([14]) is real, as can be seen from Eq. ( |i~7|) when the integral is taken along the 
imaginary uj axis. 
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FIGURES 



FIG. 1. Simple geometry for determining the effect of a thin film on a perfectly conducting plate 

FIG. 2. The numercially calculated effect of a 35 nm thick Au film on a perfectly conducting 
surface. Lower curve, no coating; middle curve, Au film; upper curve, perfectly conducting film. 
The Au film effect is of order 50% of the perfectly conducting film effect in the 100-200 nm range. 
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